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I. 



INTRODUCTION 



THE PROBLEM 

Air Traffic Control has become a safety issue of great 
importance during the last decade because of the many 
near-miss or tragic accidents that have occurred at military 
and civilian airports worldwide. The main reasons for these 
accidents are: 

* Air Traffic Control system failure, 

* Air Traffic Control erroneous procedures, 

* Pilot error, 

* Weather conditions, 

* Increased air traffic. 

This issue has even greater significance at military air 
training bases because of the very high ratio of student 
pilots to experienced pilots that are using the air space 
and because of the large volume of aircraft activity in the 
air in specific areas. It is standard military training 
practice for a large number of training aircraft to be 
assigned to the same radio channel for aircraf t-to-aircraf t 
and for aircraft-to-air traffic controller communications 
and to have a large number of training aircraft following 
approximately the same air pattern with the same air speed 
and altitude. These conditions can increase the probability 



of breakdowns in synchronization and communication between 
aircraft and between the aircraft and the air traffic 
controllers and can therefore increase the probability of 
accidents occurring. 

Kalamata Air Force Base is an Air Force training base 
located in southern Greece. This Base utilizes two types of 
aircraft for training purposes, T-37's and T-2's. The two 
different types of aircraft used at Kalamata Air Force Base 
ordinarily have different flight schedules, flight 
capabilities, and training missions that would allow them to 
take-off and begin their training flights at different 
times. However because of current operational constraints it 
is common for a group of T-2 and T-37 aircraft to complete 
their missions at approximately the same time so that they 
simultaneously return to the local Air Traffic Pattern. It 
is during the simultaneous approach of the returning 
aircraft that a critical safety problem arises due to the 
traffic congestion and pilot fatigue. Furthermore, the 
inefficient aircraft schedules interfere with the 
performance of the scheduled activities. A more efficient 
scheduling of aircraft take-offs can remedy the safety 
problems and also can diminish the necessity for aircraft to 
wait in order to get into the mission areas. 

This thesis provides a computer simulation model 
programmed in GPSS for the IBM-PC for analyzing current 
procedures for the efficient management and control of the 
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air traffic of Kalamata Air Force Base. The thesis contains 
the model logic, the GPSS program, model validation and 
model results. 
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II. MODEL DEVELOPMENT 



A. DESCRIPTION 

The environment that is modeled is the operation of the 
Kalamata Air Force Training Base in Greece. This operation 
is modeled as a sequential multiserver limited capacity 
queuing system. The server elements of the model consist of 
the base, two runways, an air traffic pattern and eight 
training mission areas. 

The base, runways and air traffic pattern are illustrated 
in Figure 2.1. The air traffic pattern is modeled as a set 
of sequential servers consisting of two entrance points (EP1 
and EP2 ) , the initial point (IP), the low initial point 
(LIP), the break point and the base key. EPl is the 
entrance point for the aircraft returning from a western 
mission area and EP2 is the entrance point for aircraft 
returning from an eastern mission area. 

The mission areas are represented in the model as single 
points. Specific mission area training activity is not 
modeled. The eight mission areas are the neighboring areas 
around the airport and are illustrated in Figure 2.2. The 
areas west of the airport are numbered one through four and 
the areas east of the airport are numbered five through 
eight . 
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Figure 2.1 Air Traffic Pattern 
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Figure 2.2 Mission Areas 
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The calling population consists of a basic training 
squadron employing T-37 type aircraft an advanced training 
squadron employing T-2 type aircraft and occasional 
aircraft from other bases. 

B. MODEL LOGIC 

The general flow through the model for a typical aircraft 
is as follows: 

A take-off time is scheduled for the aircraft from one 
of the two squadrons and the aircraft enters the base 
activity queue at that time. If there is no runway or air 
space conflict with other aircraft landing or taking-off, 
the aircraft takes-off and proceeds to the base departure 
point where it requests and is assigned to a mission area to 
carry out the scheduled training activities. 

If all training mission areas are occupied the aircraft 
is assigned to the mission training area with the smallest 
waiting line and in the case of ties the aircraft is 
assigned to the mission area with the earliest expected 
departure time for the occupying aircraft. All the aircraft 
that wait for a mission area to become available, maintain a 
maximum altitude of 7000 feet at the corresponding areas 
and perform training maneuvers consistent with the altitude 
safety requirements. Once the aircraft enters the mission 
training area it stays in the mission area for a standard 
length of time for the aircraft type and performs the 
scheduled activities. 
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At the completion of the mission area training 
activities the aircraft returns to the base and enters the 
returning air traffic pattern and follows the sequence of 
events as described below and illustrated in Figure 2.1. 

Aircraft returning from the western mission areas enter 
the air traffic pattern at entrance point EP1 and aircraft 
returning from the eastern mission areas enter the air 
traffic pattern at entrance point EP2 . The sequence of air 
legs in the base air traffic landing pattern and the 
aircraft capacity of each air leg is as follows. 

Air Leg Aircraft Capacity 



EP1-EP2 2 
EP2-IP 2 
IP-LIP 2 
LIP-Break Point 3 
Break Point-Base Key 3 
Base Key-Final 2 



If any of the first four air legs is at full capacity 
when an aircraft attempts to enter, the aircraft must orbit 
at the air leg entry point and wait for entry. If the Break 
Point-Base Key leg is at full capacity when an aircraft 
attempts to enter, the aircraft attempting entry returns to 
the LIP-Break Point air leg entry point. If the Base Key- 
Final leg is at full capacity when an aircraft attempts to 
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enter, the aircraft attempting entry performs "go around", 
that is, it cancels the landing and attempts a "close 
approach" pattern. 

Aircraft in the "close approach" pattern re-enter the air 
traffic pattern at the Base Key point if no other aircraft 
is in the Break point-Base Key air leg and if no other 
aircraft is waiting for take-off. Otherwise they re-enter 
the air traffic pattern by joining the waiting line at the 
LIP entry point. 

Upon being cleared for landing the aircraft will either 
land and return to the squadron or if flying constraints 
allow the aircraft will perform a touch and go landing. If 
a touch and go landing is made the aircraft will either 
enter the "close approach" pattern or attempt to re-enter 
the air traffic pattern with associated probabilities of 
0.20 and 0.80. If there is a waiting line at the LIP point 
or if the IP-LIP air leg is at full capacity the touch and 
go aircraft will either continue to LIP or re-enter the air 
traffic pattern at the EP1 entry point with associated 
probabilities of 0.70 and 0.30. 

It is also possible that the flow of aircraft in the air 
traffic landing pattern can be interrupted by emergency 
events. Emergency events such as engine failure, oil 
pressure failure, low fuel, or landing gear failure, are 
common occurrences that require the aircraft with the 
emergency to land as soon as possible. When such an 
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emergency event occurs, the aircraft with the emergency pre- 
empts all other aircraft in the landing pattern. During the 
emergency the nonemergency aircraft in the landing pattern 
upon reaching the IP, LIP, or Break air leg points leave the 
pattern and return to the entrance point EP1 and join the 
orbiting waiting line. If a nonemergency aircraft is on the 
Base-Final air leg when the emergency situation arises the 
aircraft continues and makes a full stop landing. 

During the emergency the aircraft present in the 
mission areas do not leave their areas but after completing 

their scheduled activities orbit at a lower altitude until 
the emergency ends to avoid interfering in the local traffic 
pattern. Aircraft returning to the air traffic pattern will 
not enter the pattern but will orbit at the entrance point 
until the emergency ends. 

If an aircraft take off is scheduled during an emergency 
event the take-off is delayed until the emergency event is 
completed. 

C. INPUT AND OUTPUT 

All input data for the model as well as suggestions 
for model structure was provided by the instructor 
pilots and from the control tower personnel at the 
Kalamata Air Force Base. The general inputs to the model 
consisted of individual aircraft characteristics and 
performance data, alternate take off schedules, mission 
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area constraints, time distributions for assignments in the 
mission areas, and the time distances between the 
reference points for the air legs in the air traffic 
pattern. 

The model outputs contain information pertaining to the 
aircraft of each squadron, the air traffic pattern and 
mission area utilization. The following output is 

available . 

1. Aircraft total flight time distribution. 

2. Number of aircraft by squadron in the east and 
the west mission areas. 

3. Number of aircraft take-offs. 

4. Number of entries in each mission area. 

5. Number of entries in each leg of the air 
traffic pattern. 

6. Maximum number of aircraft waiting in each 
queue . 

7. Average wait time for each queue. 

8. Waiting time distributions for the following 

queues : 

a. Air traffic controller. 

b. Entry points, EP1 and EP2 . 

c. Initial point, IP. 

d. Low Initial point, LIP. 
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III. GPSS APPLICATION 



A. BRIEF GPSS DESCRIPTION 

GPSS is the General Purpose Simulation System language 
developed by IBM for modeling and simulating queuing 
systems. GPSS was used to model the Kalamata Air Force Base 
air traffic operations. The GPSS program is included in 
Appendix A. 

GPSS uses the process interaction approach for modeling 
in which the model entities are either temporary or 
permanent. The temporary entities are called 
transactions and the permanent entities are called 
facilities and storages. The transactions represent the 
calling population and the facilities and storages 
represent the service centers. Transactions interact with 
other transactions and with the facilities and 
storages. In the Kalamata Air Force Base model the 
calling populations of aircraft are represented by 
transactions and the mission areas and the air traffic 
landing pattern segments are represented by facilities 
and storages. 

The modeling and programming approach in GPSS is to 
define a set of programming statements called blocks 
that represent the entrance and flow of the 
transactions into the queuing system composed of the 
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facilities and storages. There can be many 

transactions simultaneously moving through the blocks. 
At any point in time each transaction is positioned at a 
block and most blocks can hold many transactions 
simultaneously. The transfer of a transaction from one 
block to another occurs instantaneously at a specific time 
or when some change of system condition occurs. Time in 
the GPSS model is managed by the next event sequence 
with the simulation clock changing at nonuniform discrete 
time points when the state of the system changes. 
Transactions continue to move through the system until they 
either encounter a waiting line or service time delay. 

In GPSS simulated clock time is an integer value 
whose scale value is chosen by the programmer. The unit 
of time is not specifically stated but is implied by 
providing all times in the same units. In the Kalamata Air 
Force Base model the unit of time used is the second. 

B. MODEL STRUCTURE IN GPSS 

A brief description of some of the important 
programming blocks and storage areas used in the GPSS 

Kalamata Air Force Base model are contained in this 

section. 

GENERATE and TERMINATE blocks: Transactions are created 
and enter the system at one or more GENERATE blocks and 
are removed from the simulation at TERMINATE blocks. 
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The time and frequency with which transactions enter the 
system are controlled by the GENERATE block. In the 

Kalamata Air Force Base model GENERATE blocks are used for 
the entry of training aircraft from the squadrons for 
take off assignments and for the entry of occasional 
aircraft from other bases into the Kalamata air traffic 

landing pattern. 

ADVANCE block: The ADVANCE block will hold 

transactions for a specified or computed number of time 
units. The purpose of the ADVANCE block is to hold 
the transactions in service. In the Kalamata Air Force 
Base model ADVANCE blocks are used to simulate the time 
delays associated with take off delays, training mission 
areas and transit from point to point in the air leg 
segments of the air traffic landing pattern. 

TEST block: The TEST block is used to manage or transfer 
transactions based upon the test conditions. In the 

Kalamata Air Force Base model the TEST block is used to 
prevent aircraft from entering the system after the daily 
training period and to assure that the aircraft in the 

system are correctly processed in order to complete all 
landings after the daily training period ends. 

GATE block: The GATE block is used as a gate to 

interrupt the flow of transactions depending upon 
conditions that set the gate to "open" or "closed". In the 
Kalamata Air Force Base model GATE blocks are used to 
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prevent aircraft from continuing in the air traffic landing 
pattern or the take-off queue during an emergency event. 

SELECT block: The SELECT block is used to direct the 
flow of transactions. In the Kalamata Air Force Base 
model SELECT blocks are used to assign the aircraft to the 
training mission areas after take-off. 

JOIN, REMOVE, COUNT and MARK blocks: The JOIN, REMOVE, 
COUNT and MARK blocks are used to collect, remove, 
count and identify transactions in the queue. In the 
Kalamata Air Force Base model if the training mission areas 
are occupied the JOIN, REMOVE, COUNT and MARK blocks 
are used to determine current aircraft assignments based 
upon the shortest waiting lines for the areas. 

FACILITIES and STORAGE areas: GPSS FACILITIES and 
STORAGE areas are used to collect and hold transactions 
for time delays that can be associated with service or 
performance of the transactions. A GPSS FACILITY can 
hold one transaction. A GPSS STORAGE area can hold more 
then one transaction. In the Kalamata Air Force Base 
model eight FACILITIES model the eight mission areas, one 
FACILITY models emergency aircraft, one FACILITY models 
the air traffic controller, four facilities model 
aircraft synchronization, and six STORAGE AREAS have been 
used to model the air leg segments of the air traffic 
landing pattern as described in the following page. 
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FACILITY 


MODEL 


101 


Mission area 1. 


102 


Mission area 2. 


103 


Mission area 3. 


104 


Mission area 4 


105 


Mission area 5. 


106 


Mission area 6. 


107 


Mission area 7. 


108 


Mission area 8. 


DANGER 


Emergency aircraft. 


CNTR 


Traffic controller. 


DUMY1 


Aircraft synchronization 


DUMY2 


Aircraft synchronization 


DUMY3 


Aircraft synchronization 


DUMY4 


Aircraft synchronization 



STORAGE AREA 


Air Leg Segment 


CAP1 : 


EP1-EP2 


CAP2: 


EP2-IP 


CAP3: 


IP-LIP 


CAP4: 


LIP-BREAK POINT 


CAP 5: 


BREAK POINT-BASE KEY 


CAP 6: 


BASE-FINAL 
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C. GATHERING STATISTICS WITH GPSS 

GPSS automatically records data and collects queue 
statistics for transactions that pass through a storage 
area. In addition to the previously described GPSS 
programming blocks that manage the flow of transactions 
there are several block types that are specifically designed 
to gather statistics on transactions. These blocks and 
their application are described in this section. 

QUEUE and DEPART statistic blocks: QUEUE and DEPART 

blocks are used to identify specific data collection 
points in the queue. 

Data on transactions that move through the queue and 
enter and leave associated QUEUE and DEPART blocks is 

collected as separate sets of queue statistics. In the 
Kalamata Air Force Base model queue statistics were 
accumulated over the entire period of the simulation, in 
each queue of the air traffic pattern and in each queue of 
the mission areas. The accumulated statistics for each of 
the model queues are identified in the following list. 

Queue Statistics 

Maximum number of aircraft. 

Average number of aircraft 

Standard deviation of the number of aircraft. 

Average waiting time. 

Number of entries. 
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Number of entries that did not wait for entry. 

Average wait time to enter the queue. 

TABULATE, TABLE and QTABLE statistic blocks: The 

TABULATE, TABLE and QTABLE blocks are used to collect 
data for frequency and cumulative frequency tables. 

The TABLE and QTABLE blocks define the transaction 
characteristics that are to be counted and the range of 
the frequency tables. Data on characteristics of the 
transactions that pass through a TABULATE block are 
automatically collected. The TABULATE block may be used 
anywhere in the GPSS program and can collect data on all 
of the transactions computed characteristics. The 
TABULATE, TABLE and QTABLE blocks are used in the Kalamata 
Air Force Base model to tabulate the total flight time 
distribution for the two types of aircraft and waiting time 
distributions for all of the model queues. 
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IV. MODEL VALIDATION 



The Kalamata Air Force Base model was validated by 
comparing model output to historical data for a specific 
take off schedule for the squadrons of aircraft. Historical 
data was made available for each aircraft for a three day 
period. The historical data is included in Appendix C. 

For each aircraft the historical data consisted of the 
interarrival times for take off, the time spent waiting in 
the air traffic controller queue, the time spent waiting to 
enter area 3 and the total flight time. The historical data 
was accumulated and averaged for the three days of base 
activity. Using the historical interarrival rates as model 
inputs the model generated waiting times and total flight 
times were compared to the historical data. 

The historical interarrival rate appeared to be 
nonhomogeneous on a daily basis. Therefore the operations 
day was broken into time periods for which the interarrival 
rates were homogeneous. These time periods were: 



1) 


07:30 - 


08:40 


5) 


11:40 - 


14:00 


2) 


08:40 - 


09:10 


6) 


14:00 - 


14:30 


3) 


09:10 - 


11:10 


7) 


14:30 - 


16:20 


4) 


11:10 - 


11:40 


8) 


16:20 - 


17:30 
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For each time period the interarrival times appeared to 
have an exponential distribution. This hypothesis was 
tested using the Kolmogorov-Smirnov test and was not 
rejected. The hypothesis test results are included in 
appendix D. With this base information the model was run 
for a period of eight days using the hypothesized 
exponential distributions for the take off interarrival 
times. The validation comparison tests follow. 

Figure 4.1 displays a plot of the comparison of the 
actual wait time with model generated wait time for entry to 
area 3, and Figure 4.1a displays the regression line for the 
regression model : 

Simulated Data=A+B*( Actual Data). 

The slope of the regression line is .97 with a standard 
deviation of .027. The hypothesis B=1 was tested with the 
t-test and was not rejected. The probability level for the 
t-statistic was .16. 

Figure 4.2 displays a plot for the comparison of actual 
wait time with model generated wait time in the air traffic 
controller queue and Figure 4.2a displays the regression 
line for the regression model: 

Simulated Data=A+B* (Actual Data). 

The slope of the regression line is .98 with a standard 
deviation of .043. The hypothesis B=1 was tested with the 
t-test and was not rejected. The probability level for the 
t-statistic was .43. 
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Figure 4.3 displays a plot of the comparison of actual 
total flight time with model generated total flight time 
and Figure 4.3a displays the regression line for the 
regression model: 

Simulated Data=A+B*( Actual Data). 

The slope of the regression line is .96 with a standard 
deviation of .019. The hypothesis B-l was tested with the 
t-test and was not rejected. The probability level for the 
t-statistic was .09. 

These results provide a validation of the model and show 
that the model data compares favorably with the historical 
data for the specified take off schedules. 
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ACM vs SIMULATED UAITINC TIMES 



... V 
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Figure 4.1 Comparison of Actual Wait Time 
with Simulated Wait Time in Area 3 
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on Actual Wait Time in Area 3 



fiCTUBL vs SIHULftTID «IM TIKES 
(IT LINE OP 




Figure 4.2 Comparison of Actual Wait Time 
with Simulated Wait Time at Line Up 
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Figure 4.2a Regression of Simulated Wait Time 
on Actual Wait Time at Line Up 
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ML vs SIMULATED TOTAL FLIGHT TIME 
Aircraft T37 




; TIME 



Figure 4.3 Comparison of Actual Flight Time 
with Simulated Flight Time 
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Figure 4.3a Regression of Simulated Flight Time 
on Actual Flight Time 
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V. RESULTS-ANALYSIS 



This section contains a description of four basic and 
two additional take-off schedules for the squadron of 
training aircraft for Kalamata Air Force Base that were run 
in the model and a comparison analysis of the resulting air 
traffic controller congestion in order to determine a 
reasonable efficient take-off schedule. In order to run 
each take-off schedule in the model, the GPSS program had to 
be modified. These program modifications are included in 
Appendix B. 

A. SCHEDULES 

Schedule 1: This schedule consists of a forty minute 
take-off cycle. The cycle consists of two consecutive 
fifteen minute take-off periods followed by ten minutes of 
no take-off activity. The T-37 aircraft squadron assigns 
five aircraft for take-off in the first fifteen minute 
period. The T-2 aircraft squadron assigns five aircraft for 
take-off in the second fifteen minute period. This cycle is 
repeated until all the training aircraft are scheduled for 
take-off. Aircraft that can not take-off in their assigned 
period are recycled for later take-off. Aircraft taking-off 
have priority use of the runway over normally landing 
aircraft. This schedule was proposed for analysis by the 
two squadrons at Kalamata Air Force Base. 
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Schedule 1A: 



This schedule is derived from and is 



identical to Schedule 1 except that the number of aircraft 
scheduled for take-off in each fifteen minute period is four 
instead of five. This schedule was suggested as a result of 
reviewing the model output data for Schedule 1. It was 
thought that this change would decrease the number of 
training flights while also decreasing the air traffic 
congestion and waiting times. 

Schedule IB: This schedule is derived from and is 
identical to Schedule 1 with the take-off cycle period 
extended to sixty minutes. This schedule was also suggested 
as a result of reviewing the model output data for Schedule 
1. It was thought that this change would also decrease the 
number of training flights while decreasing the air traffic 
congestion and waiting times. 

Schedule 2: This schedule consists of an eight-minute 
take-off cycle. The cycle consists of two consecutive 
four-minute take-off periods. The T-37 aircraft squadron 
assigns one aircraft for take-off in the first four-minute 
period. The T-2 aircraft squadron assigns one aircraft for 
take-off in the second four-minute period. This cycle is 
repeated until all the training aircraft are scheduled for 
take-off. Aircraft that can not take-off in their assigned 
four-minute period are recycled for later take-off. This 
schedule was also suggested as a result of reviewing the 
model output data for Schedule 1. It was thought that this 
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model output data for Schedule 1. It was thought that this 
change would increase the number of training flights while 
decreasing the air traffic congestion and waiting times. 

Schedule 2A: This schedule is derived from Schedule 2 and 

is identical to Schedule 2 except that the take-off cycle is 
extended to ten minutes with two five-minute take- off 
periods. This schedule was suggested as a result of 
reviewing the model output data for Schedule 2. 

Schedule 3: This schedule is a reproduction of the 

schedule currently in use at Kalamata Air Force Base. The 
schedule contains no structure and take-off times are 
scheduled at random. The distributions of current 
interarrival times for take-off times were analyzed and 
found to be exponential for different periods during the 
day. This analysis is contained in Appendix C and was also 
used for model validation. 

B . RESULTS 

The model output data for each schedule are contained in 
Appendix E. The measures of effectiveness used in comparing 
schedules were performance and ef f iciency/saf ety . 
Performance is measured by the average number of training 
aircraft scheduled. Efficiency is measured by the average 
waiting times in the mission area queues. Safety is 
directly related to efficiency in that the smaller waiting 
times mean less hazardous flying conditions. The empirical 
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distributions for both of these measures for each schedule 



are contained in Appendix E. 

The summary results for Schedule 1 are of interest because 
this schedule was suggested by the training squadrons at 
Kalamata Air Force Base. This schedule, by comparison with 
the other schedule results, does not provide the highest 
values for performance and efficiency/safety. It was 
because of these results that the other schedules were 
derived from Schedule 1 by sensitivity analysis. 

The summary results for Schedule 3 are of interest 
because this schedule is a reproduction of the schedule 
currently in use at Kalamata Air Force Base. This schedule 
by comparison with the other schedule results, also does not 
provide the best values for performance and 
efficiency/safety. 

It was found that by applying Schedule 1A and 
comparing the results to Schedule 1 that performance 
decreased but that efficiency/safety improved drastically. 
Figure 5.1 compares the cumulative mission area waiting time 
distributions and Figure 5.2 compares the cumulative LIP 
point waiting time distributions for Schedules 1 and 1A. 

Schedule 2 results were more efficient than Schedule 1 
results. However, by applying Schedule 2A mission area 
waiting times decreased dramatically and provided the most 
efficiency/safety measure for all schedules examined. 
Figure 5.3 compares the cumulative mission area waiting time 
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distributions and figure 5.4 compares the cumulative LIP 
point waiting time distributions for schedules 2 and 2A. 
Figure 5.5 compares the cumulative mission area waiting time 
distributions for schedules 1A, 2A and 3. It is obvious 
that Schedule 2A is preferable to Schedule 1. 

Schedule 3 results were more efficient than Schedule 1 
results. The performance measure for Schedule 3 however is 
the lowest of all the schedules examined. Figure 5.6 
compares the cumulative LIP point waiting time distributions 
for schedules 1A, IB, 2A and 3. 

C. ANALYSIS 

From the above summary results the preferred schedules 
appear to be Schedules 1A and 2A. An analysis of variance 
was performed, using the function "ANOVA" from the OA3660 
APL WORKSPACE, to test the hypothesis that the mean mission 
area waiting time differences for schedules 1, 1A, 2, 2A are 
not significant. The analysis of variance results are 
contained in Table 5.1. These results show that the null 
hypothesis of no significant differences between mean 
mission area waiting times is rejected at significance 
levels greater than .995. 

The Sum of Squares from the previous analysis was broken 
into three components in order to test for individual 
effects rather than just a schedule effect using individual 
degrees of freedom. The results of this analysis are 
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contained in the tables 5.2 and 5.2a. These results show 



that the hypothesis of 


no significant difference between the 


compared mean mission 


area waiting times for the selected 


schedules is rejected 


for each comparison at significance 



levels greater than .975. 

TABLE 5.1 

ANALYSIS OF VARIANCE RESULTS 





ANOVA TABLE 


SOURCE DF 


SS MS F 


SCHEDULE 3 


232864.74 77621.59 27.77 


ERROR 4 


11181.05 2795.27 


TOTAL 7 

R-SQUARE = 0.954 


244045.78 


OVERALL MEAN = 


200.68 


TREATMENT EFFECTS 


-87. 56266. 68***. 50 14.39 


-68.5 24.02 


0.435 -17.92 


68.5 -24.02 


-0.435 17.92 
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TABLE 5.2 



ANOVA WITH INDIVIDUAL DEGREES OF FREEDOM 
SCHEDULE 1 SCHEDULE 1A SCHEDULE 2A SCHEDULE 2 



44.62 


181.62 


49.138 


443.34 


7.62 


6.75 


197.15 


232.99 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


-1 


-1 


0 


0 


0 


0 


0 


0 


0 


0 


-1 


-1 


1 


1 


1 


-1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


-1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


-1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


-1 



D = S bij 2 Z=Sb ij X j w^-Z^/Di 



8 


— 




8 


716.45 


64162 . 575 


4 


-708.48 


125485.970 


4 


415.77 


43216.173 


2 


-137.00 


9384.500 


2 


48.04 


1153.920 


2 


0.87 


0.378 


2 


1284.50 


642.250 
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TABLE 5.2a 



ANALYSIS OF VARIANCE RESULTS 
USING INDIVIDUAL DEGREES OF FREEDOM 



SCHEDULES 


SS 


DF 


MS 


F 


( 1 , 1A ) vs ( 2 , 2A ) 


64162.575 


1 


64162.575 


22.954 


1 vs 1A 


125485.970 


1 


125485.970 


44.892 


2 vs 2A 


43216.173 


1 


43216.173 


15.460 


RESIDUAL 


11181.048 


4 


2795.262 




TOTAL 


244045.760 


7 
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SCHEDULE 




(X 100) 

HIM TIME 



Figure 5.1 Cumulative Waiting Time Distribution 
in the Mission Areas 
Schedules 1 and 1A 
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SCHEIE 1 vs lfi 




BAITING TIME 



Figure 5.2 Cumulative LIP Waiting Times 
Schedule 1 and 1A 
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SCHEDULES 2 MI 2ft vs UfliTil TIME 




USIIIHC HUE 



Figure 5.3 Cumulative Waiting Time Distribution 
in the Mission Areas 
Schedules 2 and 2A 
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LIP UhITING MS 
SCHEDULE 2 vs 2A 




WAITING TIME 



Figure 5.4 Cumulative LIP Waiting Times 
Schedule 2 and 2A 
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SCHEDULES id, 2d AND 3 vs WAITING TIME 




(X 100) 

WAITING TIME 



Figure 5.5 Cumulative Waiting Time Distribution 
in the Mission Areas 
Schedules 1A 2A and 3 
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IIP UfilTING HUES COIffiEISON 
Scedules 16, 26, IS, 3 




Figure 5.6 Cumulative LIP Waiting Times 
Schedules 1A IB 2A and 3 
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VI . CONCLUSIONS - SUMMARY 



A. CONCLUSIONS 

If an added emphasis is to be placed upon the 
efficiency/safety factor in scheduling the aircraft then 
Schedule 2A is the preferred schedule. If Schedule 2A 
conflicts with local base operations due to other operating 
constraints then Schedule 1A is the next best schedule. The 
difference in mean mission area waiting times for these two 
schedules is made more significant when it is also 
realized that landing aircraft have runway priority in 
schedule 2A and do not in Schedule 1A. Figure 6.1 compares 
the cumulative mission area waiting time distributions and 
figure 6.2 compares the cumulative LIP point waiting time 
distributions for schedules 1A and 2A. 

If added emphasis is to be placed upon the performance 
factor in scheduling the aircraft then Schedule 2 is the 
preferred schedule. However, the increase of twenty four 
scheduled aircraft is at the expense of more than a ten fold 
increase in mean mission area waiting times that contributes 
to air traffic congestion and pilot fatigue and should be 
avoided if possible. Again, if Schedule 2 conflicts with 
local base operations due to other operating constraints 
then Schedule 1 is the next best schedule. It is 
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recommended that the efficiency/safety factor be the 
deciding factor in selecting schedules. 

Table 6.1 contains a summary comparison between the 
different schedules. In this table the first column 
contains the average number of flights, the second column 
contains the percentage of the aircraft that waited in the 
mission area queues, the third column contains the 
percentage of the aircraft that waited more than 180 seconds 
in the mission area queues, the fourth column contains the 
percentage of the aircraft that waited more than 30 seconds 
at LIP and the fifth column contains the average conditional 
waiting time, in seconds, in the mission area queues. 



TABLE 6.1 

SUMMARY COMPARISON BETWEEN SCHEDULES 





# of A/C 
scheduled 


% of A/c 
waiting 
in 

m. areas 


% of A/C 
waiting 
>180"in 
m. areas 


% of A/C 
waiting 
>30"in 
LIP 


Avg . Conditional 
waiting time 
in seconds 


Sch . 2A 


106 


2 


0 


0.8 


11 


Sch. 1A 


106 


7 


5.6 


1.6 


45 


Sch . 2 


130 


26 


18 


2.5 


203 


Sch. IB 


89 


7 


7 


2 


230 


Sch . 3 


99 


14 


13 


3 


368 


Sch. 1 


130 


26 


17 


3.5 


496 
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Figure 6.3 compares the cumulative mission area waiting 
time distributions and figure 6.4 compares the cumulative 
LIP point waiting time distributions for schedules 1 and 2. 



SCHEDULES IS ID 211 vs ilTIHC TIME 




UAIT1KG TIHE 



Figure 6.1 Cumulative Waiting Time Distribution 
in the Mission Areas 
Schedules 1A and 2A 
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Figure 6.2 Cumulative LIP Waiting Times 
Schedules 1A and 2A 
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SCHEDULES 1 AW 2 vs SITING TIME 




HAIM TIME 



Figure 6.3 Cumulative Waiting Time Distribution 
in the Mission Areas 
Schedules 1 and 2 
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IIP UAITINC TIMES 
Schedule i vs 2 



r 

•j 



F 

J 

0 

E 

J 

I 

S 

T 

l 




Ijl TIH5 TIHE 



Figure 6.4 Cumulative LIP Waiting Times 
Schedules 1 and 2 
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B . SUMMARY 



This thesis presents an IBM-PC GPSS model for studying 
aircraft take-off schedules at Kalamata Air Force Base in 
Greece. The model is specific to Kalamata Air Force Base 
because of the structure and sequence of the queues used in 
modeling the air traffic control points and the training 
mission areas. The GPSS program can be transferred to any 
IBM-PC compatible computer that will run GPSS. The model 
can therefore be used at the Kalamata Air Force Base to 
continue to examine other aircraft take-off schedules and 
to help reducing fuel consumption. 

Six aircraft take-off schedules were examined and a 
comparison of results was based upon factors of performance 
and efficiency/safety. The take-off schedule currently used 
at Kalamata Air Force Base is shown to be a poor performer 
with a low efficiency/safety factor. The take-off schedule 
proposed for use by the training squadron personnel at 
Kalamata Air Force Base is shown to have the highest 
performance factor and the worst efficiency/safety factor. 

The specific schedule of ten minutes of a ten minute take- 
off cycle consisting of two consecutive five minute take-off 
periods from each squadron with runway priority given to 
landing aircraft is shown to be the best schedule based upon 
the recommended emphasis placed upon efficiency/safety 
rather than on numerical performance. 
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APPENDIX A 



This appendix contains the GPSS program for the Kalamata 
Air Force Base model. 



1 0 ********************************************************************** 



20 * * 

30 * AIR-TRAFFIC CONTROL SIMULATION * 

40 * * 

50 * * 



60 *************************************************************** ******* 

70 * 



80 * 

90 * EXPONENTIAL DISTRIBUTION 

100 * 



110 EXP FUNCTION RN1,C24 

0,0/. 1 , . 104/. 2.. 222/. 3.. 355/. 4,. 509/. 5. .69/. 6. .915/. 7, 1.2/. 75. 
. 8 , 1 . 6/ . 34 , 1 . 83/ . 83 , 2. 1 1/ . 9 . 2. 3/ . 92 . 2 . 52/ . 94 , 2. 8 1 / . 95 , 2. 99/ . 9o 
.97, 3. 5/. 98, 3. 9/. 99 ,4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120 * 

1 30 * 

140 MEAN FUNCTION Cl.DS 

3600 , 380 . 1 1 /5400 .221.60/1 3200 , 488.62 

1 5000 ,242. 63/23400 , 557 . 83/25200 ,124.58 

31300,494.60/36600,221.30 

130 * 



1 

i 



j • ^ 



160 


* 








170 




GENERATE 


I” 1 

X 


•.CREATE A SINGLE TRANSACTION 


130 


SWITCH! 


LOGIC S 


: SET LOGIC SWITCh 1 


190 




ADVANCE 


300 


; T-37 TAKE OFF 


20n 




LOGIC r 


1 


; RESET LOGIC SWICH 1 


210 




ADVANCE 


2100 


; NONE T-37 TAKES GFF 


220 




TRANSFER 


,SWITCH1 




230 




GENERATE 


v 900,1 

jL 




240 


SWITCH2 


LOGIC 8 


; SET LOGIC SWICH 2 
; T-2 TAKE OFF 
•RESET LOGIC SWICH 2 


250 




ADVANCE 


300 


260 




LOGIC R 




270 




ADVANCE 


2100 


; NONE T-2 TAKES OFF 


280 




TRANSFER 


, SWITCH2 




290 




GENERATE 


90G0.FNTEXP 


; CREATE EMERGENCY EVENTS 


300 




SEIZE 


DANGE 




310 




ADVANCE 


600.180 


; EMERGENCY HOLDS 


320 




RELEASE 


DANGE 




330 

— . if* 


.4. 


TERMINATE 




; EMERGENCY TERMINATES 


35 0 




GENERATE 


9100,FNTEXP 


; A/C FROM OTHER BASES 
; ASSIGN TO A PARAMETER 


360 




ASSIGN 


i n 

1 • *- 


370 




MARK 


4' 


; MARK THE TIME 


380 




TRANSFER 


, EPI 




390 


* 








400 


* CRE 


ATE AIRCRAFTS 


FROM TWO DIFFERENT 


S0UADR3MS ON BASE 


410 


* 








420 




GENERATE 


60.10 


•.CREATE T-37 'S 

•THE SIMULATION STOPS AFTER 


430 




TEST L 


NI6NRT, 1 ,PISTA 
1 ,PISTA 


440 




GATE LS 




450 




ASSIGN 


1,1 


; ASSIGN TO A PARAMETER 


460 




TRANSFER 


,DOWN 




470 


* 








430 




GENERATE 


60.10 

NTGNRT, 1 ,PISTA 
2,PISTA 


; CREATE T-2'S 


490 




TEST L 




500 




GATE LS 
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51 0 


ASS I G\ 


1,0 


; ASSIGN "j A PARAMETER 


520 DOWN 


QUEUE 


CNTp 


•WAIT CL T O c THE RUNWAY 


Z7.C) 


GA" FV 


DANGS 


; D3NT MOVE IF EMERGENCY HOLDS 


540 * 


GATE 3E 


CAPS 


: CHECK POP A\C ON BASE LES 


LdC 


S“ I In 


CNTF. 


; CAPTURE THE CONTROLLER 


560 


DEPART 


CNTR 


•GOING FOR LINE UP 


570 LIN'UP 


AD VAN EE 


i 0 , 0 


; LINE-UP CHECK 


580 


RELEASE 


CNTR 


• TAKE-OFF 


590 


MARK 


h 


{START FLIGHT TIME r CR T-2 A/C 


600 


MARK 


7 


; START FLIGHT TIME FOR T-37 A/C 


610 


ADVANCE 


100,10 


; AFTER TAKE OFF TO DEPART POINT 


615 


TRANSFER 


. 25 , , OTHER 


; FORMAT IONS, INSTR.FL. , CPM 


620 * 








630 


SELECT E 


2,101 ,103,0, F.QUEUP 


; FIND EM D TY AREA IF EXISTS ONE 


640 AREA 


TEST E 


BV$F I RST , 1 , WEST 


; CHOOSE EAST OR WEST AREA 


650 * 








660 * E 


AST A R 


E A S 




670 * 








680 


SAVEVALUE 


13+,P1 


{ RECORDS T-37'S 


69C 


SAVEVALUE 


14+, i 


•RECORDS TOTAL # OF A/C 


700 WAIT 


QUEUE 


P2 


•WAIT IN THE QUEUE 


710 


TEST E 


PI , 1 , TTWO 


; A/C TYPE T-37 


720 


SEIZE 


P2 


; CAPTURE THE AREA 


730 


DEPART 


P2 


; LEAVE THE QUEUE 


740 


JOIN 


FARM 


; JOIN THE GROUP 


750 


MARK 


60 


; COUNT THE TIME IN THE AREA 


760 


ADVANCE 


2100,300 


•DO SCHEDULED ACTIVITY 


770 


RELEASE 


P2 


•FREE THE AREA 


780 


REMOVE 


FARM 


;OUT OF THE GROUP 


790 * 








800 


GATE FV 


DANGE 


; IF EMERGENCY ORBIT UNDER 


810 TOIP 


ADVANCE 


120,30 


; GOING TO EP1 


820 I PE 


QUEUE 


CAP! 


; ENTER THE QUEUE 


830 


SEIZE 


DUMY1 


: SYNCHRONIZE THE AIRCRAFTS 


840 


GATE FV 


DANGE 


; IF EMERGENCY ORBIT THERE 


850 


ENTER 


CAF'l 


•CAPTURE THE STORAGE 


860 


TRANSFER 


SIM, DIR, DLY 




870 DLY 


ADVANCE 


10 




880 DIR 


DEPART 


CAF'l 


; LEAVE THE QUEUE 


890 


RELEASE 


DUMY1 




900 


ADVANCE 


120, 10 


; GOING TOWARDS EP2 


910 


LEAVE 


CAP! 


•FREE THE STORAGE 


920 * 








930 BAK 


QUEUE 


CAF'2 


{ENTER THE QUEUE 


940 


SEIZE 


DUMY2 


•SYNCHRONIZE THE AIRCRAFTS 


950 


GATE FV 


DANGE 


; IF EMERGENCY ORBIT THERE 


960 


ENTER 


CAF'2 


{ENTER THE STORAGE 


970 


TRANSFER 


SIM, DRC, DLL 




930 DLL 


ADVANCE 


10 




990 DRC 


DEPART 


CAP2 


{DEPART THE QUEUE 


1000 


RELEASE 


DUMY2 




1010 


ADVANCE 


70.10 


{GOING TOWARDS IP 


1020 


LEAVE 


CAF'2 


{FREE THE STORAGE 


1030 * 








1040 RWAY 


QUEUE 


C-AP3 


{ENTER THE QUEUE AT IP 


1050 


SEIZE 


DUMY3 


{SYNCHRONIZE THE AIRCRAFTS 


1060 


GATE FV 


DANGE, EP I 


;IF EMERGENCY GO TO EPI 


1070 


ENTER 


CAP3 


{ENTER THE STORAGE 


1080 


TRANSFER 


SIM, DIREC, DELAY 




1090 DELAY 


ADVANCE 


10 




1100 DIREC 


DEPART 


CAP3 


{DEPART THE QUEUE 


1110 


RELEASE 


DUMY3 




1120 


ADVANCE 


65,5 


{GOING TOWARDS LIP 


1130 


LEAVE 


CAP3 


{FREE THE STORAGE 


1140 * 








1150 LIP 


QUEUE 


CAP4 


{ENTER THE QUEUE AT LIP 


1160 


SEIZE 


DUMY4 


{SYNCHRONIZE THE AIRCRAFTS 


1170 


GATE FV 


DANGE, EP I 


;IF EMERGENCY GO TO EPI 


1180 


ENTER 


CAP4 


{ENTER THE STORAGE 



50 



1 190 




TRAN5FEF 


SIM, EIRE, DE_A 


l^Or 


nc A 


AD7AN~- 


10 


1210 


DIRE 


DEP ART 


CA-4 


i on/ 




PE-E-ET 


D'JMY 4 






ADVANCE 


CAP4 


1240 




LEAVE 


1Z5C 




GATE F/ 


DAM-E.E=: 


1260 




GATE SNF 


CAPS, CNTU 


1270 

12S0 


* 


ENTER 


CAPS 


1290 




ADVANCE 


45.15 

CAPS 


1300 




LEAVE 


1310 




GATE FV 


DANGE.FAST 

CAP6,G0AR 


1320 


BASE 


GATE SNR 


1330 




ENTER 


CAP6 


1340 




ADVANCE 


45.15 
CAP 6 


1350 

1360 


* 


LEAVE 




1370 




GATE FV 


DANGE.TELC 


1380 




GATE FV 


cntr , Goar 


13R0 




TEST NE 


PI, 2, OUT 


1400 




-r E£ - E 


pi; vtt 
MP7, 3900, TELO 


1410 




TEST LE 


1430 


MORE 


ADVANCE 


io,! 


1433 




T ES“ NE 


PI .2, CNTU 
20+ , i 
, CNTU 


1436 




SAVEVALUE 


1440 




TRANSFER 


1450 


OUT 


TEST 5 


MP4. 900, MORE 


1460 


TELO 


ADVANCE 


20,! 


1470 




TEST NE 


PI, 2, PISTA 


1480 


FINISH 


TEST E 


Pl.l , TABUi_T2 
T0TAL37 


1490 




TABULATE 


1500 




TRANSFER 


.PISTA 


1510 


TABULT2 TABULATE 


TOT ALT 2 


1520 


PISTA 


TERMINATE 




i53:< 


* 






1540 


CNTU 


ADVANCE 


60,5 


1550 




GATE FV 


DANSE.CRIT 
.20, .CLOSE 

cap!. EPI 1 

Qf CAP4 , 0 , EP I 


1560 


FST 


TRANSFER 


1570 


HERE 


GATE SNF 


1580 


HERE1 


TEST E 


1590 




ADVANCE 


120,20 


1600 

1610 


* 


TRANSFER 


,Ll! 


1620 


EPI 


ADVANCE 


170,20 
, I Ft 


1630 

1640 


* 


TRANSFER 




1650 


GOAR 


ADVANCE 


30,10 

CAP5. CLOSE 
, HERE 


1660 




GATE SNE 


1670 




TRANSFER 


1680 


EPI1 


TRANSFER 


. 30,HERE1 ,EPI 


1690 


* 






1700 


TTWO 


SEIZE 


P2 


1710 




DEPART 


P2 


1720 




JOIN 


FARM 


1730 




MARK 


60 


1740 




ADVANCE 


2400,300 


1750 




RELEASE 


P2 


1760 




REMOVE 


FARM 


1770 

1780 


* 


TRANSFER 


,TOIP 


1790 


* 


W E 


ST AREA 


1800 


* 






1810 


WEST 


SAVEVALUE 


15+ ,P1 


1820 




SAVEVALUE 


16+ , 1 


1840 




QUEUE 


P2 


1850 




TEST E 


PI , 1 , TTWB 



; LEAVE THE QUEUE 

; GOING TO BREAK PC IN" 

; FREE THE STORAGE 

•IF EMERGENCY : GO ~C ENTR. =, C I ! ~ 

; CHECK: A/C ON DOWNLINE DON'" CLE_P 



; THEY ARE SYNCHRONIZED HERE 

; DOWNWING LEG 

;FREE THE STORAGE 

•IF EMERGENCY : GO TO ENTR.P3IN" 

; CHECK #A/C ON BASELEG DON'T OEJEU 

; ENTER THE CAPACITY 
; BASE LEG 

; FREE THE CAPACITY 

; IF EMERGENCY DO FULL STOP 
; IF A/C IN THE RUNWAY G0-AR3UNE 
; A/C FROM OTHER BASE 
; CHECK FOR A/C TYPE T-37 
; CHECK FLIGHT TIME OF T-37 
•TOUCH AND GO 

; COUNT THE NUMBER OF LANDS 

; CHECK FLIGHT TIME FOR FORE. A/C 
; FULL STOP LANDING 
; CHECK FOR OTHER TYPE A/C 
•IF A/C TYRE IS T-37 CONTINUE 
; DISTRIBUTION OF THE TOTAL TIME 

; TABULATE T-2 FLIGHT TIME 



; CONTINUE TO LIP 

; IF EMERGENCY ; GO TO EN7R . -0 1 NT 
•207. CONTINUE FOR CLOSE PATTERN 
; CHECK # A/C FROM IP TC LIP 
; IF A/C WAITS AT LIP, GO TO EPI 
; GOING TO LIP 



-.GOING TO EPI 



; GOING AROUND 

; IF NO A/C IN DOWNWING DO CLOSE 

;307. GO TO EPI 

; CAPTURE THE AREA 
.■LEAVE THE QUEUE 
: JOIN THE GROUP 
: MARK THE TIME IN THE AREA 
: DO THE SCHEDULED ACTIVITY 
;FREE THE AREA 
;OUT OF THE GROUP 



; RECORDS T-37'S 
; RECORDS TOTAL # A/C 
; ENTER THE QUEUE 

: CHECK THE TYPE O r THE A/C 



51 



:Si :> 




SEIZE 


P2 


1870 




DEBAR" 


P2 


IBS' 




J L I N 


FARM 


1 S Q 0 
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REMOVE 
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TRANSFER 
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SLCT 
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2220 




SAVEVALUE 
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2230 




REMOVE 


FARM, ALL, ,P2, X32 


2240 




TRANSFER 
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2243 
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7 TABLE 
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CNTR, 0,60, 40 


2310 


CAP1 


QTABLE 


CAP 1,0,30,40 
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101,0,120,32 


2360 


AREA2 


QTABLE 
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2370 


AREA3 


QTABLE 


103,0,120,32 


2380 


AREA4 


QTABLE 


104,0,120,32 



; A/C T-37 CAPTJRES THE AREA 
■LEAVE T HE DLEUE 
■JOIN THE GROUF 
I MARK THE TINE IN "HE AREA 
IDO THE SCHEDULED ACTIv I" t 
■FREE THE AREA 
; 0L7 CF THE GRDU= 

; GOING TO EF2 



; A/C T-2 CAPTURES THE AREA 
•LEAVE THE QUEUE 
; JOIN THE GROUP 
; MARK THE TIME IN THE AREA 
;DO SCHEDULED ACTIVITY 
•FREE THE AREA 
; OUT OF THE GROUP 
; GOING TO EP2 



; CHECK FLIGHT TIME FOR T-2 A/C 
; GO FOR FULL STOP 

; 207. GO TO I PE, OTHERS CONTINUE 
;807. GO TO I PE 

; IF NO A/C FOR TAKE OFF CONTINUE 
; DOING CLOSE PATTERN 
; BASE KEY 

; CHECK THE QUEUES IN ALL AREAS 
; CHECK IF ALL QUEUES ARE OCCUPIED 

; FIND THE EARLIEST OCCUPIED AREA 

, AREA ; REMOVE FROM THE GROUP 

; SELECT AREA WITH THE MIN QUEUE 

, AREA ; REMOVE FROM THE GROUP 

; A/C NOT GOING TO MISSION AREAS 
•PERFORM SHEDULED ACTIVITIES 
; BACK TO THE AIRPORT 



; TOTAL TIME DISTRIBUTION 
; TOTAL TIME DISTRIBUTION 
; TIME DISTRIBUTION IN CNTR QUEUE 
;TIME DISTRIBUTION IN CAP1 QUEUE 
; TIME DISTRIBUTION IN CAP2 QUEUE 
; TIME DISTRIBUTION IN CAP3 QUEUE 
; TIME DISTRIBUTION IN CAP4 QUEUE 
;TIME DISTRIBUTION IN AREA Q101 
; TIME DISTRIBUTION IN AREA C102 
;TIME DISTRIBUTION IN AREA G103 
; T I ME DISTRIBUTION IN AREA 0104 
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23V0 


ARz.42 


£T A3i_E 




± v - } 0 } a A -0 , --‘a. | i i Mz. D x 2 ; K I B_i . a Or*, i ,\ A 'An A 


2400 


AREA-6 


C~A3_E 




10t, 0,120,32 ; T I ME D I SIR I BUT I ON IN AREA 


'■> r - 


rrKiA / 


2 i AEi_h 




1 07 ,0,1 2U .32 ; TIME D I STR I BUT I ON J % ARE A 


'-t r 

*L\! 

2430 


AREAS 

* 


STABLE 




108,0, 120,32 ; TIME DISTRIBUTION IN AREA 


2440 

2450 


FIRST 

* 


BVAR I ABi_E 




( P2 ' E ' 1 0 1 +P2 ' E ' 1 02+P2 ' E ' 1 03+P2 ' E ' 1 04 ) 


2460 


* 


D E F I N 


E 


THE STORAGES 


2470 


CAP1 


ST0RA5E 




>s 


2480 


CAP2 


STORAGE 




2 


2490 


CAP3 


STORAGE 




4 - 


2500 


CAP4 


STORAGE 




T 


2510 


CAP5 


ST0RA3E 




T 


2520 

nc7 


CAP6 

* 


STORAGE 




>-i 

aL 


2540 


* 


TIME 


g 


EGKENT 


OCC 

2560 

2570 


SMRT 


GENERATE 

tcct r 

TERMINATE 




30600 ; TIME ARRIVES AT 30600 

NTDOWN , NTFINISH ;WAIT UNTIL ALL THE A/C LAND 

1 ; SHUT OFF THE RUN 
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This appendix 
for the different 



contains 
take off 



the GPSS program modifications 
schedules . 



SCHEDULE 1 



10 ********************************************************************** 
20 * * 

30 * AIR-TRAFFIC CONTROL SIMULATION * 

40 * * 

50 * * 

60 ********************************************************************** 
70 * 

SO * 

90 * EXPONENTIAL DISTRIBUTION 

100 * 

110 EXP FUNCTION RN1.C2* 

0.0/. 1 , . 104/. 2, . 222/. 3, . 355/. 4 , . 509/ 

.6,1.6/. 84,1. 8*/. S3, 



•- % t ‘ ^ ! « ~T « * 

12/. 9,2.3/. 92, 



69/. 6,. 915/. 7. 1.2/. 75, 1.3E 



’7u2/ . 94 



l. 81/. 95,2. 99/. 96,3. 2 



u- | i ■ U / VJ •_* / • CU I i. • X j— • / i i- ■ / • / X i X ■ UX. / ■ / ~ ^ X i U i / a /U a Xi 

97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995. 5. 3/. 998, 6. 2/. 999, 7/. 9996, 8 
120 * 

130 * 

140 MEAN FUNCTION Cl.DB 



3600, 



380.11/5400,22 

IT A A 



1 . 60/13200,483. 61 



1 5000 , 242 . 63/23400 , 557 . 83/ 25260 , 1 24 . 58 
3 1 800 , 494 . 60/36600 ,221. 30 
150 * 



170 




GENERATE 




; CREATE A SINGLE TRANSACT I OF 


180 


SWITCH1 


LOGIC S 


; SET LOGIC SWITCH 1 


190 




ADVANCE 


300 


•T-37 TAKE OFF 


200 




LOGIC R 


1 

X 


; RESET LOGIC SWICH 1 


210 




ADVANCE 


2100 


•NONE T-37 TAKES OFF 


220 




TRANSFER 


,SU1ITCH1 


230 




GENERATE 


^,900, 1 




240 


SWITCH2 


LOGIC S 


: SET LOGIC SWICH 2 


250 




ADVANCE 


300 


:T-2 TAKE OFF 


260 




LOGIC F: 


»-k 

X 


; RESET LOGIC SWICH 2 


270 




ADVANCE 


2100 


; NONE T-2 TAKES OFF 


280 




TRANSFER 


.SWITCH2 
9000, FNf EXP 




290 




GENERATE 


; CREATE EMERGENCY EVENTS 


300 




SEIZE 


DANGE 




310 




ADVANCE 


600.180 

DANSE 


; EMERGENCY HOLDS 


320 




RELEASE 




330 




TERMINATE 




; EMERGENCY TERMINATES 


340 


* 








350 




GENERATE 


9100,FN$EXP 


: A/C FROM OTHER BASES 


360 




ASSIGN 


1 •> 
1 , a. 


: ASSIGN TO A PARAMETER 


370 




MARK 


4 


; MARK THE TIME 


380 




TRANSFER 


,EPI 




390 


* 








400 


* CREATE AIRCRAFTS 


FROM TWO DIFFERENT 


SQUADROMS ON BASE 


410 


* 








420 




GENERATE 


60,10 


: CREATE T-37'S 


430 




TEST L 


N$GNRT , 1 ,PISTA 


; THE SIMULATION STOPS AFTER 


440 




GATE LS 


1 ,PISTA 
1,1 




450 




ASSIGN 


; ASSIGN TO A PARAMETER 


460 




TRANSFER 


,D0WN 




470 


* 






480 




GENERATE 


60.10 


: CREATE T-2'S 



SCHEDULE 1A 



10 ****** ********* ********************* ***** *************** ************** 



20 * + 

30 * AIR-TRAFFIC CONTROL SIMULATION * 

40 * * 

50 * * 



60 ********************************************************************** 



70 * 

80 * 

90 * EXPONENTIAL DISTRIBUTION 

100 * 

110 EXP FUNCTION RN1,C24 
0,0/. 1 , . 104/ . 2, . 222/. 3 , . 355/ . 4 , . 509/ . 5, . 69/. 6 , . 915/. 7 . 1 . 2/. 75, 1 . 3 
.8, 1.6/. 84, 1.83/. 88, 2. 12/. 9,2. 3/. 92,2. 52/. 94,2. 81/. 95,2. 99/ . 96,3. 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120 * 

130 * 



-•o 

•y 



140 MEAN FUNCTION C1,D8 
3600 . 380 . 1 1 /540Q ,221. 60/ 1 3200 , 488 . 62 
1 5000 , 242 . 63/23400 , 557 . 83/25200 ,124.58 
3 1800,494. 60 /3660C , 22 1 . 30 

150 * 



160 * 



170 




GENERATE 


} ,,i 


; CREATE A SINGLE TRANSAC" 


180 


SWITCH1 


LOGIC S 


; SET LOGIC SWITCH 1 


190 




ADVANCE 


300 


; T-37 TAKE OFF 


200 




LOGIC P 


i 


•RESET LOGIC SWICH 1 


210 




ADVANCE 


2100 


; NONE T-37 TAKES OFF 


220 

230 




TRANSFER 

GENERATE 


,SWITCH1 
, ,900,1 


240 


SWITCH2 


LOGIC S 


A 7 

jL 


; SET LOGIC SWICH 2 


250 




ADVANCE 


300 


; T-2 TAKE OFF 
•RESET LOGIC SWICH 2 


260 




LOGIC R 


9 


270 




ADVANCE 


2100 


; NONE T-2 TAKES OFF 


280 




TRANSFER 


■SWITCH2 
9000 ,FN$EXP 
DANBE 


290 

300 




GENERATE 

SEIZE 


; CREATE EMERGENCY EVENTS 


310 




ADVANCE 


600,180 


; EMERGENCY HOLDS 
; EMERGENCY TERMINATES 


320 

330 

340 


* 


RELEASE 

TERMINATE 


DANBE 



350 
360 
370 
380 
390 * 
400 * 



GENERATE 


9100, FNf EXP 


; A/C FROM OTHER BASES 


ASSIGN 


1,2 


; ASSIGN TO A PARAMETER 


MARK 


4 


; MARK THE TIME 


TRANSFER 


, EPI 





CREATE AIRCRAFTS FROM TWO DIFFERENT SQUADROMS ON BASE 



f!K 



410 * 

420 

430 

440 

450 



GENERATE 75.10 

TEST L N$GNRT, 1 ,PISTA 

GATE LS 1 ,PISTA 

ASSIGN 1,1 



; CREATE T-37'S 

;THE SIMULATION STOPS AFTER 

; ASSIGN TO A PARAMETER 



460 
470 * 
480 
490 



TRANSFER , DOWN 

GENERATE 75,10 : CREATE T-2'S 

TEST L N$GNRT, 1 ,PISTA 



500 

510 

520 DOWN 
530 
540 * 



GATE LS 
ASSIGN 
QUEUE 
GATE FV 
GATE SE 



2,PISTA 

1,0 

CNTR 

DANSE 

CAP6 



•.ASSIGN TO A PARAMETER 
; WAIT OUT OF THE RUNWAY 
■DONT MOVE IF EMERGENCY HOLDS 
•CHECK FOR A\C ON BASE LEG 



JJ 



SCH=.DLi_E 1 E 



10 ********************************************************#**#********-#* 



20 * * 

30 * AIR-TRAFFIC CONTROL SIMULATION * 

40 * * 

50 * * 



60 ********************************************************************** 
70 * 

80 * 

90 * EXPONENTIAL DISTRIBUTION 

100 * 

110 EXP FUNCTION RN1,C24 

0 , 0/ . 1 , . 1 04/. 2 , . 222/ . 3 , . 355/ . 4 . 509/ . 5 , . 69/ . 6 , . 91 5/ . 7 , 1 . 2/ . 75 , 1 . 36 
.8, 1.6/. 84, 1.83/. 88, 2. 12/. 9,2. 3/. 92,2. 52/. 94,2. 81/. 95,2. 99/. 96,3. 2 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120 * 

130 * 

140 MEAN FUNCTION C1,D8 

3600 , 380 .11/5400,221. 60/1 3200 , 43S . 62 
1 5000 , 242 . 63/23400 , 557 . 83/25200 ,124.58 
3 1 800 , 494 . 60/36600 ,221. 30 
150 * 

1 60 * 



170 




GENERATE 


? ’ ’ 1 
1 


; CREATE A SINGLE TRANSACTION 


180 


SWITCR1 


LOGIC S 


; SET LOGIC SWITCH 1 


190 




ADVANCE 


300 


; T-37 TAKE OFF 


200 




LOGIC R 


1 


; RESET LOGIC SWICH 1 


210 




ADVANCE 


3300 


; NONE T-37 TAKES OFF 


220 




TRANSFER 


, SWITCH1 


230 




GENERATE 


*,900,1 




240 


SWITCH2 


LOGIC S 


; SET LOGIC SWICH 2 


250 




ADVANCE 


300 


;T-2 TAKE OFF 


260 




LOGIC R 


2 


; RESET LOGIC SWICH 2 


270 




ADVANCE 


3300 


; NONE T-2 TAKES OFF 


280 




TRANSFER 


.SWITCH2 

9000,FNfEXP 




290 




GENERATE 


; CREATE EMERGENCY EVENTS 


300 




SEIZE 


DANGE 




310 




ADVANCE 


600, 180 
DANSE 


; EMERGENCY HOLDS 


320 




RELEASE 




37.0 




TERMINATE 




; EMERGENCY TERMINATES 


340 


* 








350 




GENERATE 


9100,FN$EXP 


; A/C FROM OTHER BASES 


360 




ASSIGN 


1,2 


■ASSIGN TO A PARAMETER 


370 




MARK 


4 


; MARK THE TIME 


o*8 0 




TRANSFER 


,EPI 




390 


* 








400 


* CREATE AIRCRAFTS 


FROM TWO DIFFERENT 


SQUADRONS ON BASE 


410 


* 








420 




GENERATE 


60.10 

N$GNRT,1,PISTA 
1 ,PISTA 


; CREATE T-37'S 


430 




TEST L 


; THE SIMULATION STOPS AFTER 


440 




GATE LS 




450 




ASSIGN 


1,1 


; ASSIGN TO A PARAMETER 


460 




TRANSFER 


, DOWN 




470 


* 








480 




GENERATE 


60.10 


; CREATE T-2'S 


490 




TEST L 


N$GNRT, 1 ,PISTA 




500 




GATE LS 


2, PI ST A 




510 




ASSIGN 


1,0 


; ASSIGN TO A PARAMETER 


520 


DOWN 


QUEUE 


CNTR 


; WAIT OUT OF THE RUNWAY 


530 




GATE FV 


DANGE 


: DONT MOVE IF EMERGENCY HOLDS 


540 


* 


GATE SE 


CAP6 


; CHECK FOR A\C ON BASE LEG 
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SCHEDULE 



10 **************************************************************-#-* +•*** *-* 



20 * + 

30 * AIR-TRAFFIC CONTROL SIMULATION * 

40 * * 



50 * * 

60 *********************************************************************** 
70 * 

80 * 

90 * EXPONENTIAL DISTRIBUTION 

100 * 

110 EXP FUNCTION RN 1,024 

0.0/. 1,. 104/. 2.. 222/. 3,. 355/. 4,. 509/. 5,. 69/. 6,. 915/. 7, 1.2/. 75, 1.38 
.8. 1.6/. 84. 1.83/. 88,2. 12/. 9,2. 3/. 92,2. 52/. 94,2. 81 /. 95 .2. 99/. 96,3. 2 
. 97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, B 
120 * 



130 * 

140 MEAN FUNCTION Cl.DB 
3600 , 380. 11/5400,221. 60/13200 , 488. 62 
15000,242. 63/23400,557. 83/25200 , 124. 58 
3 1 800 , 494 . 60/36600 ,221. 30 
150 * 

160 * 

290 GENERATE 

300 SEIZE 

310 ADVANCE 

320 RELEASE 

330 TERMINATE 

340 * 



9000, FNf EXP 
DANGER 
600.180 
DANGER 



350 
360 
370 
380 
390 * 



GENERATE 9100,FN*EXP 

ASSIGN 1,2 

MARK 4 

TRANSFER ,T0EP1 



; CREATE EMERGENCY EVENTS 

; EMERGENCY HOLDS 

; EMERGENCY TERMINATES 

; A/C FROM OTHER BASES 
; ASSIGN TO A PARAMETER 
•MARK THE TIME 



400 

410 

411 

412 

413 

414 

415 

416 

417 
420 

430 

431 
450 
460 
470 
480 

490 

491 
510 
520 
530 
540 
550 

560 

561 

562 
570 
580 
590 
600 



* CREATE AIRCRAFTS FROM TWO DIFFERENT SQUADROMS ON BASE 



* 





GENERATE 


120,10, ,4 




ASSIGN 


i , i 




TRANSFER 


, DOWN 




GENERATE 


120,10,60,4 




ASSIGN 


1,0 




TRANSFER 


, DOWN 


r 


GENERATE 


480,60,1500 




TEST L 


NfGNRT , 1 ,PISTA 




ASSIGN 


1,1 


JL 


TRANSFER 


, DOWN 


A 


GENERATE 


480,60,1740 


JL 


TEST L 


N$GNRT,1,PISTA 




ASSIGN 


1,0 


DOWN 


QUEUE 


CNTR 




GATE FV 


DANGER 




GATE SE 


CAP6 




SEIZE 


CNTR 


JL 


DEPART 


CNTR 


* 


l : 


I N E UP 


LINEUP 


ADVANCE 


70,20 




RELEASE 


CNTR 




MARK 


6 




MARK 


7 



; CREATE THE FIRST 4 T-37'S 
•ASSIGN TO A PARAMETER 

; CREATE THE FIRST 4 T-2'S 
; ASSIGN TO A PARAMETER 



; CREATE T-37'S 
; THE SIMULATION STOPS AFTER 
; ALL THE A/C HAVE LANDED 
; ASSIGN TO A PARAMETER 



CREATE T— 2 ' S 

THE SIMULATION STOPS AFTER 
ALL THE A/C HAVE LANDED 
ASSIGN TO A PARAMETER 
WAIT OUT OF THE RUNWAY 
DONT MOVE IF EMERGENCY HOLDS 
CHECK FOR A\C ON BASE LEG 
CAPTURE THE CONTROLLER 
GOING FOR LINE UP 



; LINE-UP CHECK 
; TAKE-OFF 

; START FLIGHT TIMEFOP T-2 A/C 
; START FLIGHT TIME FOP T-37 £ 



/C 
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SCHEDULE 2A 



10 **************^************************************************* ****** 



20 * * 

30 * AIR-TRAFFIC CONTROL SIMULATION * 

40 * * 

50 * * 



60 ********************************************************************** 
70 * 

QO * 

90 * EXPONENTIAL DISTRIBUTION 

100 * 

110 EXP FUNCTION RN1,C24 

0,0/. 1,. 104/. 2,. 222/. 3, .355/. 4, .509/. 5,. 69/. 6,. 915/. 7, 1.2/. 75, 1.38 
.8, 1.6/. 84, 1.83/. 88, 2. 12/ . 9,2. 3/. 92,2. 52/. 94,2. 81/ . 95,2. 99/. 96,3. 2 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120 * 

130 * 

140 MEAN FUNCTION C1,D8 

3600 , 380 . 1 1 /5400 , 221 . 60/13200 , 488. 62 
1 5000 , 242 . 63/23400 , 557 . 83/25200 ,124.58 
3 1 eOO , 494 . 60/36600 ,221. 30 
150 * 

160 * 



290 

300 




GENERATE 

SEIZE 


9000, FNf EXP 
DANGER 


; CREATE EMERGENCY EVENTS 


310 

320 

330 

340 


* 


ADVANCE 

RELEASE 

TERMINATE 


600,180 

DANGER 


; EMERGENCY HOLDS 
; EMERGENCY TERMINATES 


350 




GENERATE 


9 100, FNf EXP 


; A/C FROM OTHER BASES 


360 




ASSIGN 


1,2 


; ASSIGN TO A PARAMETER 


370 

380 

390 


* 


MARK 

TRANSFER 


4 

,T0EP1 


; MARK THE TIME 


400 

410 


* CREATE AIRCRAFTS 

* 


FROM TWO DIFFERENT 


SQUADROMS ON BASE 


411 




GENERATE 


120,10, ,4 


; CREATE THE FIRST 4 T-37'S 


412 

413 




ASSIGN 

TRANSFER 


11 

,6own 

120,10,60,4 


; ASSIGN TO A PARAMETER 


414 




GENERATE 


; CREATE THE FIRST 4 T-2'S 


415 




ASSIGN 


lf>OWN 


; ASSIGN TO A PARAMETER 


416 

417 


* 


TRANSFER 


420 




GENERATE 


600,60,1500 


; CREATE T-37'S 


430 

431 


* 


TEST L 


NfGNRT , 1 ,PISTA 


; THE SIMULATION STOPS AFTER 
; ALL THE A/C HAVE LANDED 
; ASSIGN TO A PARAMETER 


450 

460 

470 


* 


ASSIGN 

TRANSFER 


IftoWN 


480 




GENERATE 


600,60,1800 


; CREATE T-2'S 


490 

491 


* 


TEST L 


NfGNRT, 1 ,PISTA 


; THE SIMULATION STOPS AFTER 
; ALL THE A/C HAVE LANDED 


510 




ASSIGN 


1,0 


: ASSIGN TO A PARAMETER 


520 


DOWN 


QUEUE 


CNTR 


; WAIT OUT OF THE RUNWAY 


530 




GATE FV 


DANGER 


: DONT MOVE IF EMERGENCY HOLDS 


540 




GATE SE 


CAP6 


: CHECK FOR A\C ON BASE LEG 


550 




SEIZE 


CNTR 


; CAPTURE THE CONTROLLER 


560 

561 


* 


DEPART 


CNTR 


; GOING FOR LINE UP 


562 


* 


L I 


N E UP 




570 


LINEUP 


ADVANCE 


70,20 

CNTR 


; LINE-UP CHECK 


580 




RELEASE 


; TAKE-OFF 


590 




MARK 


6 


•START FLIGHT TIMEFOR T-2 A/C 


600 




MARK 


7 


; START FLIGHT TIME FOR T-37 A/C 
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20 * * 

30 * AIR-TRAFFIC CONTROL SIMULATION * 

40 * * 

50 * * 

60 ********************************************************************** 
70 * 

80 * 

90 * EXPONENTIAL DISTRIBUTION 

100 * 

110 EXP FUNCTION RN 1 , C24 

0,0/. l f . 104/. 2,. 222/. 3 ..355/. 4, .509/. 5,. 69/. 6.. 915/. 7, 1.2/. 75,1. 38 
.8, 1.6/. 84, 1.83/. 88, 2. 12/ . 9,2.3/ . 92,2.52/ .94,2.81/. 95.2.99/ .96,3.2 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120 * 

130 * 

140 MEAN FUNCTION C1,D8 

3600 , 380 . 1 1 /5400 ,221.60/1 3200 ,488.62 
15000 , 242.63/23400 , 557 . 83/25200 ,124.58 
31800 , 494 . 60/36600 ,221. 30 
150 * 



160 


* 








290 




GENERATE 


9000.FNTEXP 

DANGER 


; CREATE EMERGENCY EVENTS 


300 




SEIZE 




310 




ADVANCE 


600,180 


; EMERGENCY HOLDS 


320 




RELEASE 


DANGER 




330 




TERMINATE 




; EMERGENCY TERMINATES 


340 


* 






350 




GENERATE 


9100,FN$EXP 


; A/C FROM OTHER BASES 


360 




ASSIGN 


1,2 


; ASSIGN TO A PARAMETER 


370 




MARK 


4 


; MARK THE TIME 


380 




TRANSFER 


,T0EP1 




390 


* 








400 


* CREATE AIRCRAFTS 


FROM TWO DIFFERENT 


SQUADROMS ON BASE 


410 


* 








420 




GENERATE 


FNf MEAN , FNfEXP 
NfGNRT.l.PISTA 
.49, ,ACFTT2 


; CREATE AIRCRAFT 

•THE SIMULATION STOPS AFTER 


430 




TEST L 


435 




TRANSFER 




450 




ASSIGN 


,6own 


; ASSIGN TO A PARAMETER 


460 




TRANSFER 




470 


* 








510 


ACFTT2 


ASSIGN 


1,0 


; ASSIGN TO A PARAMETER 


520 


DOWN 


QUEUE 


CNTR 


; WAIT OUT OF THE RUNWAY 


530 




GATE FV 


DANGE 


;DONT MOVE IF EMERGENCY HOLDS 


540 




GATE SE 


CAP6 


: CHECK FOR A\C ON BASE LEG 


550 




SEIZE 


CNTR 


; CAPTURE THE CONTROLLER 


560 




DEPART 


CNTR 


; GOING FOR LINE UP 


570 


LINUP 


ADVANCE 


70.20 

CNTR 


; LINE-UP CHECK 


580 




RELEASE 


; TAKE-OFF 


590 




MARK 


6 


; START FLIGHT TIMEFOR T-2 A/C 


600 




MARK 


7 


; START FLIGHT TIME FOR T-37 A/C 


610 




ADVANCE 


100,10 


; AFTER TAKE OFF TO DEPART POINT 


615 




TRANSFER 


.25,, OTHER 


; FORMAT IONS, INSTR.FL. , CPM 


620 


* 








630 




SELECT E 


2. 101 . 108.0.F.QUEUP : FIND EMPTY AREA IF EXISTS ONE 


640 


AREA 


TEST E 


BV$F I RST , 1 , WEST 


; CHOOSE EAST OR WEST AREA 


650 


* 






660 


* E 


AST A R 


E A S 




670 


* 








680 




SAVEVALUE 


13+,P1 


; RECORDS T-37'S 


690 




SAVEVALUE 


14+ , 1 


; RECORDS TOTAL # OF A/C 


700 


WAIT 


QUEUE 


P2 


; WAIT IN THE QUEUE 
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APPENDIX C 



This appendix contains the historical data for a three day 
period. For each aircraft the historical data consisted of the 
interarrival times for take off, the time spent waiting in th< 
air traffic controller queue, the time spent waiting to entea 
area 3 and the total flight time of T-37 aircraft. All times ar< 
stated in seconds. 



Interarrival times for take off 



INARRDAYl ! Interarrival Tinea for take off Day 1) 

557 213 254 399 968 48 16 91 211 490 484 622 43 11 7 541 41 171 160 177 191 435 
472 82 27 646 305 1410 154 490 196 479 358 1637 207 84 568 548 1211 511 
31 225 85 565 41 149 778 279 157 362 575 1147 305 288 475 383 635 118 844 
274 1623 117 173 6 31 58 24 57 133 775 92 7 142 260 7 481 77 260 587 77 
39 378 225 235 786 70 163 200 79 322 22 59 118 298 70 3 280 146 



INARRDAY2 (Interarrival Tines for take off Day 2) 

291 228 259 597 89 70 351 1415 147 179 1006 388 167 43 478 22 135 18 21 357 257 
606 207 413 17 260 95? 414 182 1424 486 33 213 304 1188 362 261 177 211 
91 53 310 14 128 267 304 90 37 274 11 804 959 580 425 1491 280 887 6 235 
947 114 155 141 157 46 63 152 155 151 113 128 37 1101 1142 919 242 132 
1460 216 397 61 1029 393 26 221 80 52 534 85 99 261 382 60 112 



IKARKBAY3 ! Interarrival Tines for take off Day 3) 

438 170 1353 577 124 144 08 338 260 636 330 27 547 15 427 52 127 42 267 23 11 
?4 323 141 596 115 798 421 148 171 657 214 20 371 106 1356 78 120 22 51 
114 64 ?22 1135 481 1832 665 296 45 34 519 230 204 14 102 708 1238 1134 
1027 270 5t4 916 221 49 333 7 168 87 91 2 504 252 35 62 39 55 144 165 67 
107 2 116 i46 357 270 980 1422 869 334 1097 873 277 447 312 621 274 30 21 
492 71 452 530 741 143 676 
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Time spent waiting in the air traffic controller queue 



Waiting Controller queue Area 3 

Time dayl day2 day3 dayl day2 day3 



0 


42 


45 


49 


14 


12 


13 


60 


42 


38 


31 


0 


0 


1 


120 


12 


20 


11 


1 


1 


0 


180 


3 


5 


4 


0 


1 


0 


240 


2 


2 


0 


0 


0 


1 


300 


0 


0 


0 


0 


1 


0 


420 


0 


2 


0 


0 


0 


0 


480 


0 


0 


0 


0 


0 


1 


540 


0 


0 


1 


0 


0 


0 


660 


0 


0 


0 


1 


0 


0 


720 


0 


0 


0 


1 


0 


0 


840 


0 


0 


0 


0 


0 


1 


960 


0 


0 


0 


1 


0 


0 


1020 


0 


0 


0 


0 


1 


0 



Total flight time aircraft T-37 



Time 


dayl 


day2 


day3 


3960 


6 


5 


8 


4080 


23 


18 


13 


4200 


12 


23 


17 


4320 


10 


9 


8 


4440 


4 


3 


2 


4560 


0 


2 


1 


4800 


0 


1 


0 
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APPENDIX D 



This appendix contains the summary statistics and the 
hypothesis test results for the eight time periods for which 
the interarrival times for take off were homogeneous. These 
times periods were: 



1) 


07:30 


- 08:40 


5) 


11:40 - 


14:00 


2) 


08:40 


- 09:10 


6) 


14:00 - 


14:30 


3) 


09:10 


- 11:10 


7) 


14:30 - 


16:20 


4) 


11:10 


- 11:40 


8) 


16:20 - 


17:30 



The Tables D.l and D.2 contain the summary statistics 
for each time period and the figures D.l through D.8 contain 
the distribution fitting and the Kolmogorov-Smirnov test 
results . 
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Table D.l Summary Statistics 



Variable: INTERVAL! INTERVALS 


Sample size 
Average 
tied i an 
Mode 

Geo me trie me an 

Variance 1 

S land a r d cl e v i a t i on 

Standard error 

Min imam 

Max i mum 

Range 

Lower quarti 1 e 
U p p e r q u artile 
1 n t e r q u. a r 1 1 1 e r a n q e 
S Vir-wness 

Standard i zed skewness 
K u r t o s i s 

S t a n d a rdi z e d k u rta s i s 


35 

380 .114 

260 

259 

248.624 

20438 

347.041 
58 . 6607 
16 
1415 
1399 
144 
547 
403 

1 . 65105 
3.98766 
2 . 56701 
3 . 09996 


45 

221 .689 
160 
11 

116.221 
44135.7 
210 . 085 
31.3176 
7 

798 

791 

42 

413 

371 

0.928257 

2.54214 

-0.118849 

-0.162741 


V a i i a b .1 e : 


INTERVALS 


I NTERVAL4 


Cample s.i ze 


47 


27 


A zerage 


488.617 


242 . 63 


Median 


305 


204 


Mode 


304 


14 


Geometric m e a n 


284.774 


144.641 


Variance 


229164 


49612.2 


Standard d o v .i. a 1 i a r i 


478.711 


222 . 738 


S tandard e r r o r 


69.8271 


42.866 


Minimum 


20 


14 


Mar. imum 


1832 


778 


Range 


1812 


764 


Lower quarti. 1 e 


120 


45 


Upper quarti le 


646 


304 


Iri ter quarti le range 


526 


259 


Skewness 


1 . 28888 


1.1451 


Standard i zed skewness 


3 . 60734 


2.42911 


f urtosis 


0.696499 


0.434659 


Standardized kur tosis 


0.974683 


0.461025 



63 



Table D.2 Summary Statistics 



Variabl e : 


INTERVALS 


INTERVAL 6 


S a mi p 1 p s i r. e 


36 


39 


A V P r a Q P 


557.833 


124.897 


Med ian 


450 


107 


hade? 


383 


p 


Geometric mean 


322.289 


69.1587 


Variance 


190050 


20045 . 8 


S t and a r d d e v i a t i on 


435.947 


141 . 583 


Stan d ard e r ro r 


72.6579 


22.6714 


Minimum 


6 




Maximum 


1623 


775 


Range 


1617 


773 


L o w e r q u a r t i .1 e 


228 


39 


Upper quartile 


901 . 5 


155 


T n terquarti 1 e range 


673 . 5 


116 


Skewness 


0 . 705837 


3 . 12258 


S t a n d ard i eel s k ewness 


1.72894 


7.96105 


Kur tosis 


-0 . 301500 


12.1983 


S t a n ( J a r d i z e d k u r t o s i s 


-0 . 36927 


15.5498 


Variable : 


I NTERVAL7 


INTERVALS 


Samp 1 e si ze 


38 


30 


Average 


494.605 


221.3 


Med ian 


345.5 


130.5 


Mode 


77 


112 


G e o m e t r i c m e a n 


327.136 


129.011 


Variance 


165600 


42845 


Standard deviation 


406 . 94 


206 . 99 


S Land ard e r r o r 


66.0143 


37.7911 


M j n .7. mum 


26 


3 


M-i x imum 


1460 


741 


Range 


1434 


738 


Lowe r q ua r tile 


216 


70 


Upper quart! 1 e 


869 




In terquarti le range 


653 


r? S3 


SI ewness 


0.90161 


1 . 10583 


S t an d ardize d skewness 


2 . 269 


2.4727 


Ku rtosis 


-0.289484 


0.246429 


Standard i zed kur t os is 


-0.364259 


0.275516 
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Frequency Histogram 




Estimated KOLMOGOROV statistic 
Estimated KOLMOGOROV statistic 
Estimated overall statistic DN 
Approximate sign 1 f icance 1 eve 1 



DPLUS 
DMINUS 
= 0.09 
- 1 



0 . 0733633 
0.0923791 
791 



Figure D.l Distribution Fitting Interval 1 
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Frequency Histogram 




Estimated KOLMOGOROV statistic 
Estimated KOLMOGOROV statistic 
Estimated overall statistic DN 
Approximate significance level 



DF'LUS = 0.13479 7 
DMINUS = 0.111456 
= 0.134797 
= 0.386891 



Figure D.2 Distribution Fitting Interval 2 
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0.4 0,8 



< i 

llfc 



1,6 



Estimated 1' OL MOSQROV statistic 
Estimated KOLMOGOROV statistic 
Estimated overall statistic DN 
Approximate significance level 



DF'LUS = 0 . 0708765 
DMINUG =• 0.0722221 
= 0.0722221 
= 1 



Figure D.3 Distribution Fitting Interval 3 
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P Vi "•VI 

nr" 



U\ r I •, nr, •, ft*. 

ain 




Ee t i ma ted KOLMOGOROV statistic 
Estimated KOLMOGOROV statistic 
Estimated overall statistic DN 
A p p r o x j in a t e s j. g n i f i c; a n c e 1 eve 1 



DF’LLJS - 0.0935019 
DM I NUB = 0.0871474 
= 0.0935019 
= 1 



Figure D.4 Distribution Fitting Interval 4 
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i IM C 



Estimated KOLMOGOROV statistic DPLUS = 0-0638073 
Estimated KOLMOGOROV statistic DMINUS = 0.110322 
Estimated overall statistic DN = 0.110322 
Approximate significance level = 0.999951 



Figure D.5 Distribution Fitting Interval 5 
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Estimated KOLMOGOROV statistic 
Estimated KOLMOGOROV statistic 
Estimated overall statistic DM 
Approximate significance level 



DPLUS = 0.128466 
DMINIJS = 0.0882622 
= 0.128466 
- 0.540454 



Figure D.6 Distribution Fitting Interval 6 
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Estimated 
Estimated 
Es t 1 mated 



KOLMOGOROV statistic 
KOLMOGOROV statistic 
overall statistic DN 



Approxi mate significance level 



DPLUS = 0 - 0533982 
DMINUS = 0.117 
= 0.117 
= 0.999736 



Figure D.7 Distribution Fitting Interval 7 
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Frequency Histogram 




Estimated KOLMOGOROV statistic 
Estimated KOLMOGOROV statistic 
Estimated overall statistic DN 
Approximate significance level 



DF'LUS = 0.0867169 
DM1 MUG = 0.0760762 
= 0.0867169 
= 1 



Figure D.8 Distribution Fitting Interval 8 
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APPENDIX E 



This appendix contains the GPSS program outputs for 
each take off schedule. 

SCHEDULE 1 



FACILITY 


ENTRIES 


; UTIL. 


AVE. TIKE 


AVAILABLE 


OWNER 


PEND 


INTER 


RETRY 


EELAY 


101 


14 


0.B77 


2130. B6 


1 


0 


0 


0 


0 


0 


102 


14 


0. 650 


2113.79 


1 


0 


0 


0 


0 


0 


103 


j. *- 


0.747 


2166.25 


1 


0 


0 


0 


0 


0 


104 


12 


0.831 


2412.67 


1 


0 


o 


0 


0 


r j 


105 


11 


0.694 


2195.62 


1 


0 


0 


6 


0 


6 


106 


i n 


0 . BOB 


2345.33 


1 


0 


0 


0 


0 


0 


107 


12 


0.799 


2317.33 


1 


0 


0 


o 


0 


0 


103 


9 


0.594 


2300.00 


1 


0 


0 


6 


6 


0 


DANGE 


'-i 


0.041 


726.00 


1 


o 


0 


0 


0 


0 


CNTR 


130 


0.266 


71.36 


1 


o 


0 


o 


6 


f> 


DU V Y1 


116 


0.011 


3. 35 


1 


0 


0 


o 


0 


f) 


DUMY2 


160 


0.000 


0. 10 


1 

A 


6 


6 


0 


0 


6 


DUMY3 


160 


0.000 


0.00 


1 


0 


0 


0 


0 


o 


DUMY4 


624 


0.047 


2.67 


1 


0 


0 


0 


0 


0 



QUEUE 


MAX 


CCNT. 


ENTRIES 


ENTRIES (0) 


AVE.CONT. 


AVE. TIME 


AVE. (-0) 


RETRY 


101 


1 


0 


14 


11 


0. 10 


255.07 


1190.33 


0 


102 


1 


o 


14 


13 


0. 03 


71.43 


1000.00 


0 


1 03 


1 


0 


12 


11 


0.01 


29.92 


359 . 00 


0 


104 


1 


0 


4 n 
lii 


7 


0.04 


10- . 50 


262.80 


0 


105 


1 


6 


11 


9 


0.02 


53.00 


291.50 


0 


106 


1 


o 


12 


7 


0.01 


41.92 


100. 60 


0 


107 


1 


0 


12 


7 


0.05 


<* crcr <->cr 

-A. W — > ■ A. J 


372.60 


0 


103 


1 


0 


9 


5 


0.04 


157.44 


354 . 25 


0 


CNTR 


n 

A- 


0 


130 


TT 


0.09 


24.61 


32.93 


6 


CAF’l 


2 


0 


116 


109 


0.01 


3. 62 


63.29 


o 


CA D 2 


1 


0 


160 


159 


0.00 


0. 10 


16.00 


0 


CAP3 


1 


0 


160 


160 


0.00 


0. 00 


0.00 


0 


CAF’4 


"7 


0 


624 


557 


0.06 


3. 11 


26.94 


0 



STORAGE 


CAP. 


REMAIN. 


MIN. 


MAX. 


ENTRIES 


AVL. 


AVE. C. 


UTIL. 


RETRV 


DELAY 


CAP1 


2 


2 


0 


2 


116 


1 


0.40 


0.201 


0 


0 


CAP2 


2 


n 


0 


o 

z. 


160 


1 


0.32 


0. 161 


0 


0 


CAP3 


2 


2 


0 


n 

j- 


160 


1 


0.30 


0. 150 


0 


0 


CAP4 


—r 

O 


~T 

o 


0 


3 


624 


1 


1.28 


0.425 


0 


0 


CAP5 


3 


3 


0 


3 


624 


1 


0.81 


0.269 


0 


0 


CAP6 


2 


2 


0 


2 


648 


1 


0.84 


0.418 


0 


0 
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"''ABLE 


MTA*.’ 

. ’t 


c-- -jry 


RETRY RAN3E 




E;mi irk 


EUK. /. 


r-kIT 

Lfv : n 


*1 i. J - 


2 _ i ' 
















— 


0 




25.38 








0 - 


60 


§7 


92.31 








60 - 


120 


1 u 


100. oc 


CAPl 


3.e2 


17.01 


0 














- 


0 


109 


93.97 








30 - 


60 


4 


97.41 








60 - 


90 


1 


98.28 








90 - 


120 


1 


99. 14 








120 - 


150 


1 


100.00 


CAP2 


0. 10 


1.26 


0 


0 


159 


99.37 








0 - 


30 


1 


100.00 


CAF'3 


0 . 00 


0 . 00 


0 














— 


0 


160 


100.00 


CAPA 


3. 1 1 


9.92 


0 














- 


0 


557 


89.26 








0 - 


30 


38 


95.35 








30 - 


60 


29 


100. 00 


TDTAL37 


4159.52 


164.86 


0 

3B40 - 


4080 


26 


39.39 








4080 - 


4320 


31 


86.36 








4320 - 


4560 


6 


95.45 








4560 - 


4800 


3 


100.00 


T0TAi_T2 


4750.58 


133.22 


0 

4320 - 


4560 


4 


6.25 








4560 - 


4800 


36 


62.50 








4800 - 


5040 


21 


95.31 








5040 - 


5280 


T 


100.00 



TABLE 


MEAN 


STD. DEV. 


RETRY 


RANGE 




FREQUENCY 


CUM. V. 


AREA1 


255.07 


510.60 


0 




0 


i * 

X X 


78.57 








960 


- 


1080 


1 


85.71 








1200 


- 


1320 


1 

X 


92.86 








1320 


— 


1440 


1 


100.00 


AREA2 


71.43 


267.26 


0 




0 


13 


92.86 








960 


— 


1080 


1 


100.00 


AREA3 


29.92 


103.63 


0 




0 


11 


91.67 








240 


- 


360 


1 


100.00 


AREA4 


109.50 


181.31 


0 




0 


7 


58.33 








0 


- 


120 


2 


75.00 








240 


- 


360 


1 


83.33 








360 


— 


480 


2 


100.00 


AREA5 


53.00 


138.98 


0 




0 


9 


81.82 








120 


- 


240 


1 


90.91 








360 


- 


480 


i 


100. OC 
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rASLE 


tfEAN 


STD. DEV, 


RETRY 


RANGE 




r EEC2EV2 v 


~ xl 'y 7 


AREAS 


^1.92 


63. 32 


0 
















0 


- 


120 


“T 


83 ! 33 








120 


— 


24 : 


r> 


4 ~\ * '* 
X - * ■ V 


AREA7 


155. 25 


-y-r 
C • / 


0 




0 


7 


58.33 








0 


- 


120 


1 


66 . 67 








240 


- 


360 


n 


33.33 








360 


- 


480 


1 


91.67 








600 


— 


720 


1 


100.00 


AREA8 


157.44 


283.67 


0 




0 


5 


55.56 








0 


- 


120 


X 


66 . 6 7 








120 


- 


240 


2 


88.89 








840 


- 


960 


1 


100.00 



T 



1A 






FACILITY 


PMTD rrc 


JT T 1 


AVE. TIME As 


/AILAE^E 


OWNER 


PcNi.- 


iN'ThS 


FvL R\ 




101 


1 3 


0.775 


2066.31 


1 


0 


o 


0 


0 


0 


102 


13 


0.792 


2111.69 


1 


0 


0 


0 


f) 


6 


103 


11 


0.740 


2330. IB 


1 


0 


0 


0 


0 


0 


104 


11 


0.735 


2316.09 


1 


0 


0 


0 


0 


0 


105 


9 


0 . 600 


2311.22 


1 


0 


0 


0 


0 


0 


106 


9 


0. 607 


2338.22 


1 


0 


0 


0 


6 


0 


107 


6 


0.417 


2410.33 


1 


0 


0 


0 


0 


c 


10S 


6 


0.395 


22S1.83 


i 

X 


0 


0 


0 


0 

0 


0 


DANSE 


1 


0.022 


772.00 


1 


0 


0 


0 


o 


CNTR 


106 


0.20B 


67.97 


1 


0 


0 


0 


0 


0 


DUMY1 


93 


0.011 


4.33 


1 


0 


0 


0 


6 


0 


DUMY2 


123 


0 . 000 


0. 11 


1 


0 


0 


0 


0 


C) 


DUMY3 


123 


0.000 


0. 12 


1 


0 


0 


0 


0 


6 


DUMY4 




0.019 


1.28 


1 


0 


0 


o 


0 


0 



QUEUE 


MAX 


CQNT. 


ENTRIES 


ENTRIES (0) AVE.C0NT. 


AVE. TIME 


AVE. (-05 


RETRY 


101 


1 


0 


13 


13 


0.00 


0.00 


0.00 


0 


102 


1 


0 


13 


13 


0.00 


0.00 


0.00 


6 


103 


•t 

± 


0 


11 


10 


0.00 


3.45 


38.00 


f) 


104 


1 


0 


11 


10 


0.00 


0. 18 


2. 00 


0 


105 


1 

X 


0 


9 


9 


0.00 


0.00 


0.00 


6 


106 


1 


0 


9 


8 


0.01 


26.56 


239 . 00 


o 


107 


* 

X 


0 


6 


6 


0.00 


0.00 


0.00 


0 


1 08 


« 


0 


h 


r 

. i 


0.00 


13.00 


78.00 


0 


CNTR 


1 

X 


0 


106 


78 


0.01 


2.72 


10.29 


0 


CAP! 


o 


0 


93 


87 


0.01 


4.71 


73.00 


0 


CAF'2 


* 

X 


0 


4 *~>~r 


4 

1 Z.Z. 


0.00 


0.11 


1 4 . 00 


n 


CAP3 


1 


0 


123 


4 

1 JLjL 


0 . 00 


0. 12 


15.00 


6 


CAP4 


n 


6 


522 


494 


0.02 


1.33 


24.82 


o 



STORAGE 


CAP. 


REMAIN. 


MIN. 


MAX. 


ENTRIES 


AVL. 


AVE. C. 


UTIL. 


RETRY 


DELAY 


CAP1 


2 




0 


<-» 

jL 


93 


1 


0.32 


0. 162 


0 


0 


CAF'2 


2 


*"> 

X. 


0 


2 


123 


1 


0.25 


0. 123 


o 


6 


CAP3 


n 


o 


o 




123 


1 


0.23 


0.116 


0 


o 


CAP 4 




*T 


0 


7 


522 


1 


1.06 


0.355 


0 


0 


CAPS 


3 


~T 


0 


3 


522 


1 


0.67 


0.225 


0 


6 


CAF'6 


2 


2 


0 


2 


572 


1 


0.74 


0.372 


6 


r. 



TABLE 


MEAN 


STD. DEV. 


RETRY RANBE 




FREQUENCY 


CUM. 7. 


CNTR 


2.72 


5.60 


0 














— 


0 


78 


73.58 








0 - 


60 


28 


100.00 


CAP1 


4.71 


18.51 


0 














— 


0 


87 


93.55 








30 - 


60 


2 


95.70 








60 - 


90 




oq qq 








90 - 


120 


i 


1 00 ! 00 



76 



i 



£?!_£ 


MEAN 


S T D. DEV. 


RET D Y 


RA\"3E 




rr REC , JE-3 v 


r” ;►.* ■/ 


CA c 2 


c. :: 


1 . 26 


o 




o 




7 c . 19 








0 


— 


30 




1 00 . 00 


C-AF'3 


0. 12 


TC 

A . 


0 




f 




CD ' c 








0 


- 


30 


^ T 


100. 00 


CAP4 


1.33 


6.27 


0 




0 


4R4 


94. 6 £ 








0 


- 


30 


19 


98.28 








30 


— 


60 


9 


100.00 


T0TAL37 


4128.00 


136.31 


0 
















3840 


- 


4080 


25 


48.08 








4080 


- 


47.20 


nn 


90.39 








4320 


~ 


4560 


5 


100. 00 


T0TALT2 


4703.52 


128.07 


0 

4320 




4560 


6 


11.1 1 








4560 


- 


4800 


37 


79.63 








4800 


— 


5040 


9 


96.30 








5040 


— 


5280 


o 

jL. 


100.00 



TABLE 


MEAN 


STD. DEV. 


RETRY RANGE 




FREQUENCY 


CUM . 7. 


AREA1 


0. 00 


0.00 


0 


0 


13 


100.00 


AREA2 


0. 00 


0.00 


0 














— 


0 


13 


100.00 


AREA3 


3.45 


11.46 


0 














- 


0 


10 


90.91 








C) - 


120 


1 


100. 00 


AREA4 


0. IS 


0.60 
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